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ABSTRACT 


A study of Han River fish assemblage dynamics for 4 years was conducted. From April 2005 to August, 2008, 
fishes inhabiting two sites of Han River were sampled for identification. For further analysis, 40 individuals 
of the dominant species were sampled monthly from March 2006 to November 2008. The fish assemblage at 
site 1 was dominated by Zacco platypus (32.69%), while the subdominant species were Acheilognathus yamat- 
sutae (14.4%), Acanthorhodeus gracilis (9.43%), Squalidus japonicus coreanus (6.84%), and Tridentiger bre- 
vispinis (5.18%). The most abundant species at site 2 was Korean Chub (Zacco koreanus) with relative abund- 
ance of 62.45% and followed by Pungtungia herzi (10.29%), Coreoperca herzi (8.67%), and Coreoleuciscus 
splendidus (6.82%) as the subdominant species. At both sites, the endemics populations show an increasing 
pattern during the whole survey period, while the natives were declining in the last two years. 
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INTRODUCTION 


Han River with the length of 470 km, located in the geograph- 
ic center of the Korean Peninsula, is a primary water resource 
for drinking, irrigation, navigation, and recreation to more 
than 23 million people. It has 35 tributaries with two major 
branches, the Namhan River and the Bukhan River (Ministry 
of Environment Republic of Korea, 2002). 

Due to the rapid development in the economy and indus- 
tries around the central part of Korea since 1960s, the water 
quality of the Han River gradually worsened (Ministry of 
Environment Republic of Korea, 2002). Although there has 
been some improvement as a result of a great investment in 
the management of the mainstream water quality since the 
late 1980s, there are still serious pollution problems (Kim et 
al., 2007). It means that continuous monitoring of Han River 
ecosystem health is urgent and important. Since the fish is 
known to be the best indicator of such purpose, regular mon- 
itoring effort should be carried out. Short term studies on fish 
diversity in Han River has been conducted by several resear- 
chers, but work for longer period in several years has not 
been presented. 


MATERIALS AND METHODS 


Study area 


The locations of the sampling sites are presented in Fig.1. 
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Site 1 

It is located in Sinjang-dong, Hanam-si, Gyeonggi-do, 
below Paldang Bridge with the coordinates of E 127° 09’ 
26” N 37° 34’05” 


Site 2 
It is located in Deode-ri, Gapyeong-gun, Gyeonggi-do 
with the coordinates of E 37° 56’48” S 127° 28’59” 


Sampling and identification 


All fishes inhabit both sites were sampled twice a year for 
identification from 2005 until 2008. A further analysis was 
conducted to the dominant species based on the annual data. 
Forty individuals of Zacco platypus and Zacco koreanus 
were sampled monthly from March 2006 until November 
2008, except for the freezing months. 

Fish samples for community analysis were collected using 
throwing nets (mesh size 7 X 7mm), hand nets (mesh size 4 
x4mm), and gill nets (mesh size 15 x 15 mm), while for 
dominant species analysis, the fish were collected using throw- 
ing nets (mesh size 7 X 7mm). At the sampling sites, collect- 
ed fish were directly fixed in 10% formalin solution and trans- 
ported to the laboratory. Further process was done at the 
laboratory where the fish were identified using the methods 
described by Kim and Park (2002). The classification system 
followed Nelson (1994). 


Community analysis 


The dominance index were calculated using the Mc Naughton 
(1967)’s dominance index, while the Diversity index and 
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Evenness index calculations were done using the Shannon- 
Wiener function (Pielou, 1966). 
The condition factor was calculated based on analysis 


Fig. 1. Map of study area showing locations of the sites in Gye- 
onggi-do. 


described by Anderson and Neumann (1996): K=W/L? x 
10°” here K as the condition factor, W as the total weight 
and L as the total length. GSI was calculated as: GSI (%)= 
gonad weight/body weight x 100 (Strange, 1996). 


Statistical analysis 


The descriptive statistics, F test and Post-Hoc test of One- 
way ANOVA was used to check any significant dissimilarity 
of the physico-chemical parameters. 

Hierarchical Cluster Analysis was used for classifying the 
species variation found in both sites. To check the correlation 
among communityanalysis indices (Richness index, Diversity 
index, Dominance index, and Evenness index) with any of 
the physico-chemical parameters of the waters at each site, 
the Pearson’s correlation test was used. 

Linear Regression analysis was used to analyze the Length- 
weight relationship and condition factor of the dominant 
species for each site, while Avedev means were used to deter- 
mine the spawning period of the dominant species of each 
site. 
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Fig. 2. Physico-chemical Analysis of each site from Jan 2005 to Dec 2008. 
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RESULTS AND DISCUSSION 


Physico-chemical parameters 


The physico-chemical data presented in Fig. 2 shows that at 
Site 1, the range of temperature, pH, DO, BOD, SS, TN, TP 
levels are 12.9-14, 7.8-8.4, 10.5-11, 1.1-1.4, 7.1-10.2, 1.9- 
2.7, and 0.04-0.05, respectively, while at Site 2, 11.7-12.8, 
7.2-7.8, 10.3-11.9, 0.5-0.7, 0.5-0.9, 1.86-1.97 and 0.03-0.04. 
The One way ANOVA test result shows that the variation 
of temperature and DO for 4 years period at Site 2 was higher 
than at Site 1, while the variation of pH, BOD, SS, TN, and 
TP at Site 2 was lower than at Site 1. The results also sug- 
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gest that pH, BOD, SS, TN and TP at Site 1 and Site 2 are 
significantly different (P< 0.001) By comparing the means 
shown by descriptive statistics result, the levels of pH, BOD, 
SS, TN, and TP at Site | were significantly higher than at Site 
2. In addition, the F value confirms that the highest difference 
between Site 1 and Site 2 was at the levels of BOD and SS 
(F value: 71.36, 96.93, also P values respectively). 

According to EPA (1991), all the water quality parameters 
of both sites were within the range of healthy ecosystem 
standards, except for the TN level which was > 1 mg/L. 

The Site 1 F test result during survey period presents 
significant difference only at pH and TN levels (Sig. levels: 
0.037; 0.008, respectively). However, the Post Hoc test 


Table 1. Fish Assemblage and its Relative Abundance (Accumulated data of 2005-2008) 


St. 1 St. 2 

Species RA (%) Species RA (%) 
Cyprinidae 87.33 Cyprinidae 86.12 
Cyprinus carpio® 0.65 Pungtungia herzi° 10.29 
Carassius auratus° 2.16 Pseudopungtungia tenuicorpa** 4.34 
Carassius cuvieri® 0.07 Coreoleuciscus splendidus* 6.82 
Rhodeus uyekii* 0.65 Hemibarbus mylodon** 0.62 
Rhodeus notatus° 3.02 Gobiobotia brevibarba** 0.62 
Acheilognathus yamatsutae* 14.4 Microphysogobio longidorsalis* 0.74 
Acheilognathus rhombeus° 0.07 Zacco koreanus* 62.45 
Acheilognathus gracilis* 9.43 Zacco platypus° 0.25 
Pseudorasbora parva’? 2.66 Cobitidae 2.48 
Pungtungia herzi° 1.08 Misgurnus anguillicaudatus° 0.12 
Sarcocheilichthys variegatus wakiyae* 0.22 Koreocobitis rotundicaudata* 2.11 
Sarcocheilichthys nigripinnis morii* 0.14 ksookimia koreensis* 0.25 
Squalidus gracilis majimae* 0.29 Siluridae 0.50 
Squalidus japonicus coreanus* 6.84 Silurus microdorsalis* 0.50 
Hemibarbus labeo° 4.46 Amblycipitidae 1.86 
Hemibarbus longirostris° 0.5 Liobagrus andersoni* 1.86 
Pseudogobio esocinus°® 2.59 Centropomidae 8.67 
Abbottina rivulariso° 0.14 Coreoperca herzi* 8.67 
Microphysogobio yaluensis* 0.29 Odontobutidae 0.37 
Microphysogobio jeoni* 1.44 Odontobutis platycephala* 0.37 
Rhynchocypris oxycephalus° 0.5 
Zacco platypus° 32.69 
Opsariichthys uncirostris amurensis° 2.95 
Hemiculter leucisculus° 0.07 
Cobitidae 0.14 
Misgurnus anguillicaudatus° 0.07 
Koreocobitis rotundicaudata* 0.07 
Centropomidae 0.07 
Siniperca scherzeri° 0.07 
Centrarchidae 0.14 
Micropterus salmoides® 0.14 
Odontobutidae 0.65 
Odontobutis interrupta” 0.65 
Gobiidae 11.66 
Chaenogobius urotaenius°® 0.29 
Rhinogobius giurinus° 2.09 
Rhinogobius brunneus° 4.1 
Tridentiger brevispinis° 5.18 
#: Endemics, °: Natives, ©: Introduced, RA: Relative Abundance 
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Fig. 3. RA of Native, Endemics, and Introduced Species at both sites during the surveys. 


result for Site 1 only confirms that TN levels in 2008 were 
significantly lower compared to the TN levels in 2005 (Sig. 
level: 0.009) and the TN levels in 2007 (Sig. level: 0.023). 

As for Site 2, the significant variations during 2005-2008 
were at the level of pH and TP with P value of 0.013 and 
0.041, respectively. The Post Hoc test confirms difference 
for pH level between 2006 and 2008 (Sig. level of 0.01) with 
negative mean difference, while the difference for TP levels 
of 2007 and 2008 was significant (Sig. level of 0.013) with 
positive mean difference. 


The fish assemblage 


Of the 1389 total individuals, 33 species from 6 families 
were recorded during the survey period at Site 1. As shown 
in Table 1, the fish assemblage was dominated by Zacco 
platypus (32.69%), while the subdominant species (with RA 
>5%) were Acheilognathus yamatsutae (14.4%), Acanthor- 
hodeus gracilis (9.43%), Squalidus japonicus coreanus (6.84 
%), and Tridentiger brevispinis (5.18%). The annual data also 
presents Pale Chub (Zacco platypus) as the most abundant 
species with relative abundance of 49.04% (2006), 55.29% 
(2007), and 32.89% (2008), respectively. Although the 2008 
data shows that the most abundant was Acheilognathus yamat- 
sutae (23.3%), Zacco platypus with relative abundance (RA) 
of 21.2% was also considered as the dominant species since 
the RA gap between them was too small. The results showing 
that Pale Chub (Zacco platypus) as the dominant species at 
Site 1 was expected since this species was also the most 
abundant species in most of large river systems of South 
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Korea (Jang et al., 2002b) and the subdominant at Mt. Jiri 
and Mt. Seorak National Parks (Jang et al., 2002a). This 
species was also included as an abundant species in Miho 
Stream (Son, 1983), Tamjin River (Hwang and Choi, 1995) 
and Hwangryong River System (Song and Park, 1994). 

It is even considered as a representative species across the 
whole of East Asia (Chen, 1982; Shen et al., 1993). 

From Site 2, 15 species from 6 families were recorded 
during 2005-2008 (Table 1). The most abundant species was 
Korean Chub (Zacco koreanus) with relative abundance of 
62.45% and followed by Pungtungia herzi (10.29%), Core- 
operca herzi (8.67%), and Coreoleuciscus splendidus (6.82%) 
as the subdominant species. The annual data during the sur- 
veys shows that Korean Chub was always present as the dom- 
inant species. This species was distinguished from the Dark 
Chub (Zacco temmincki) after 2005 by Kim et al. (2005). 

Hemibarbus mylodon, one of the Natural Monuments, 
was surprisingly found at this area with a relative abundance 
of 0.62 in 2 years of survey (2005 and 2008). In addition, 
two endangered species (Pseudopungtungia tenuicorpa and 
Gobiobotia brevibarba) were also found at this site with 
relative abundance of 4.34 and 0.62, respectively. 

This founding was also previously reported by Jang et al. 
(2008) in the same tributary of Han River (North Han River). 
Coates et al. (2003) estimated that globally there are 155, 
154, and 372 critically endangered, and vulnerable freshwater 
fish species, respectively. Therefore, the founding of these 3 
species is important. 

The fish assemblage distribution status (native, endemics, 
introduced) during the survey period for both sites is shown 
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Dendrogram using average linkage (Between groups) 
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Rescaled distance cluster combine 


Case 0 5 10 15 20 25 
Lable Num.*777777777 eer mae Be Tarasa EE ae j 

R. giurinus 31 
R. brunneus 32 
C. urotaenius 30 
T. brevispinis 33 
O. interrupta 29 
P. esocinus 5 
H. labeo 11 
H. longirostris 12 
O. amurensis 16 
P. parva 6 
P. herzi 

C. auratus 

C. carpio 10 
A. yamatsutae 2 
S. coreanus 23 
A. gracilis 4 
Z. platypus 24 
C. cuvieri 9 
M. salmoides 28 
S. morii 21 
S. majimae 22 
M. yaluensis 15 
S. wakiyae 20 
K. rotundicaudata 25 
M. jeoni 14 
R. uyekii 17 
A. rhombeus 3 
H. leucisculus 13 
A. rivularis 1 
M. anguillicaudat 26 
S. scherzeri 27 
R. notatus 18 


R. oxycephalus 19 


Fig. 4. Similarity Hierarchical Clustering of Fish assemblage at site 1 using Complete Linkage (between groups), based on the level 
of abundance, the presence frequency in years of survey, distribution status, and the family. 


in Fig. 3. At Site 1, the graph shows that the native fish 
were the most abundant with RA of 51.42% while the 
endemics were the next abundant species (34.41%). In 
contrast, the endemics (RA of 89.34%) at Site 2 were 
strongly dominated the natives (RA of 10.66%). At both 
sites, the endemics show an increasing pattern during the 
whole survey period, while the natives were declining in the 
last two years of survey. 

The introduced species or the exotic species, each year 
contributed small number of relative abundance at Site 1 
(0.64%, 0%, 1.34%, 0.81%, respectively) and none (RA of 
0%) at Site 2. 

The highest abundance of native species was reported (Jang 
et al., 2002a, b) and confirmed at the present study. The 
endemic species (the second most abundant) in this study 
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was also reported by De Silva et al. (2007) for South Korean 
water systems, as part of the total 49 endemic species in South 
Korea (Kim, 1995). According to Jang et al. (2003), the Kore- 
an endemic species of freshwater fish are mostly distributed 
in headwater/mountain streams. This statement is relevant 
to the current results showing that the abundance level of 
endemicity at the site located in mountain area was higher 
than the one located in the lower area. 

Fig. 4 shows the hierarchical clustering for the fish assem- 
blage based on the level of abundance, the presence frequency 
in years of survey, the fish distribution status and the family 
they belong to. The dendrogram identifies five major groups 
labeled as Group A, B, C, D, and E (Fig. 4). Group A repre- 
sents the fish species which are natives/endemics, abundant/ 
less abundant/rare, present in two/three/four years of survey 
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Dendrogram using average linkage (Between groups) 


Rescaled distance cluster combine 
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Fig. 5. Similarity Hierarchical Clustering of Fish assemblage at site 2 using Complete Linkage (between groups), based on the level 
of abundance, the presence frequency in years of survey, distribution status, and the family. 
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Fig. 6. Community Analysis Indices for each site during the years 
of survey (2005-2008). 


and belong to Gobiidae family. Fish which are natives/exotic, 
less abundant, present in three/four years of survey, and 
belong to Cyprinidae are represented in Group B. Group C 
represents endemics/natives fish, most abundant/abundant, 
present in two/three/four years of survey, and belong to 
Cyprinidae. Group D represents exotic fish, rare, and only 
found in a year of survey, and belong to Cyprinidae/Centrar- 
chidae, while endemic/native fish, rare/less abundant, present 
in one/two years of survey, and belong to either Cyprinidae/ 
Cobitidae/Centropomidae are represented in Group E. 

At Site 2, the cluster presents 4 major groups labeled as 
Group A, B, C, and D (Fig. 5). Group A represents the fish 
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Fig. 7. Community Analysis Indices during the years of surveys 
(2005-2008). 


Korean J. Syst. Zool. 26(2), 105-114 


Fish Assemblage Dynamics in the Han River 


Table 2. Pearson Correlation test for per-survey Physico-chemical Parameters and Community Analysis Indices 


Richness Diversity Dominance Evenness 
St. 1 St. 2 St. 1 St. 2 St. 1 St. 2 St. 1 St. 2 

Pearson Cor 0.067 0.759* —0.414 0.463 0.588 —0.379 —0.818* —0.107 
TN Sig. (2-tailed) 0.876 0.029 0.308 0.248 0.125 0.354 0.013 0.801 

N 8.000 8.000 8.000 8.000 8.000 8.000 8.000 8.000 

Pearson Cor 0.315 = —0.318 = 0.554 — —0.869** = 
TP Sig. (2-tailed) 0.448 = 0.443 = 0.154 E 0.005 = 

N 8.000 = 8.000 = 8.000 7 8.000 = 


**Correlation is significant at the 0.01 level (2-tailed). *Correlation is significant at the 0.05 level (2-tailed). 


species which are endemics, less abundant/abundant/most 
abundant, present in either two/three/four years of survey, 
and belong to Cyprinidae/Cobitidae/Centropomidae. Natives 
and abundant fish, present in all years of survey, and belong 
to Cyprinidae are represented in Group B. Group C represents 
endemic fish, rare/less abundant, present in either one/two/ 
three years of survey, and belong to Cobitidae/Odontobutidae/ 
Siluridae/Amblycipitidae. The last group (Group D) repre- 
sents native fish, rarely found, present in one/two years of 
survey, and belong to Cyprinidae/Cobitidae. 


Community analysis 


The accumulative data of Community Analysis Indices for 
each site in 2005-2008 is presented in Fig. 6. Site 1 has much 
higher Diversity Index compared to Site 2 (1.95, 1.17, respec- 
tively), while Site 2 shows a higher Dominance Index than 
Site 1 (0.76, 0.55, respectively). 

The changes of community indices from year of 2005 to 
2008 are presented in Fig. 7 (Site 1 and 2). At Site 1, the diver- 
sity indices in 2005-2008 were 1.87, 1.47, 2.33, 2.14, respec- 
tively, indicating that in the last 2 years of survey (2007-2008) 
they were relatively increased. On the contrary, at Site 2 (Fig. 
7) the diversity indices show a significant decreasing pattern 
(1.61, 1.32, 0.94, 0.92, respectively) while the dominance 
indices were slightly increased (0.61, 0.73, 0.86, 0.79, res- 
pectively). 

The Pearson Correlation test was conducted to check any 
correlation between the physico-chemical parameters and 
the community analysis indices. The results shown in Table 
2 and 3 are only the ones with significant correlations. 

At Site 1 (Table 2), the evenness index was significantly 
correlated to TN and TP (Sig. value of 0.013 and 0.005, res- 
pectively). Both of the relationships show negative correla- 
tions, which mean the TN and TP levels in the water body 
were decreased along with the increasing Evenness indices. 
At Site 2 (Table 2), the results show a negative significant 
correlation between Richness indices and the levels of TN 
(with Sig. value of 0.029) which could indirectly affect the 
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Fig. 8. Pale Chub LWR (log transformation) and Korean Chub 
LWR (log transformation) during the period of survey (2005-2008). 


fish diversity. 

The Physicochemical data shows signs of excessive TN 
levels which leading to eutrophication although the DO levels 
did not suggest this. The results showing that the annual 
mean of TP levels at Site 1 were all above 0.025 mg, also 
indicate that there was a tendency of eutrophication (Danilov 
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Table 3. Linear Regression Statistics of LWR (Log Transformation) of Pale Chub And Korean Chub during the survey period 


Coefficients (a) 


Unstandardized 


Standardized 


coefficients coefficients t Sig. 
Model B Std. error Beta 
Pale Korean Pale Korean Pale Korean Pale Korean Pale Korean 
Chub Chub Chub Chub Chub Chub Chub Chub Chub Chub 
1 (Constant) -5.616 -5.916 0.022 0.025 -253.057 -236.94 .000 .000 
Log TL 3.287 3.4417 0.012 0.014 0.993 0.992 271.468 251.421 .000 .000 


a. Dependent Variable: Log W 


0.00 0.20 0.40 0.60 0.80 1.00 
K (Condition factor) 


0.00 0.20 0.40 0.60 0.80 1.00 
K (Condition factor) 


Fig. 9. Monthly Mean of Pale Chub and Korean Chub Condition 
Factor. 


and Ekelund, 2000). Therefore, the declining of TP may have 
been the cause of the increasing evenness indices and given 
a contribution in the increase of fish diversity at the site. 

It is possible that the uniformity in the distribution of indi- 
viduals among species were increasing due to the decreasing 
of the exceeded nutrient levels. The significant negative 
correlation between the TN levels and richness indices 
could indirectly affect the fish diversity at Site 2. However, 


112 


there is not much information available explaining the respon- 
se of fish communities in large rivers to the cascade of effects 
caused by an imbalance in nutrients and any indirect and pro- 
bably unpredictable effect of enrichment (Miltner and Rankin, 
1998). 


Dominant species population 


The Log-transformed Linear of the Pale Chub is presented 
in Fig. 8, which shows a positive allometric growth with re- 
gression coefficient b of 3.287 and Correlation coefficient 
r’=0.986. And as presented in Fig. 8, the Log-transformed 
Linear of the Korean Chub LWR also shows a positive allo- 
metric growth with regression coefficient b of 3.441 and 
Correlation coefficient r°=0.983. The statistical analysis for 
the Linear Regression of both species is presented in Table 
3 showing the significance level (P value) of < 0.001. 

A value significantly larger or smaller than 3.0 indicates 
allometric growth. A value less than 3.0 shows that that the 
fish becomes lighter (negative allometric) or greater than 3.0 
indicates that the fish becomes heavier (positive allometric) 
for a particular length as it increases in size (Wooton, 1998). 
Wooton (1998)’s statements clearly suggest that within the 
Pale Chub and Korean Chub population collected in the pre- 
sent study, the allometric coefficients indicate both species 
have positive allometric growth. 

Accumulated monthly data shows that the pale chub best 
condition was in May, June, and July which is in late spring 
to early summer (Fig. 9), while the Korean Chub best condi- 
tion was in May to July (late spring to summer) (Fig. 9). 

Monthly variations of the Pale Chub males and females 
Gonadosomatic Index (GSI) showed peaks in May to July 
which indicates the spawning season of this species at the 
particular site (Fig. 10). This result is similar with the found- 
ing by Nakamura (1952) which stated that Pale Chub spawn- 
ing period occurs from May to June. 

Fig. 10 presents the GSI Monthly variations of the Korean 
Chub both for males and females. The graph shows peaks in 
May and June indicating that the species spawns in this par- 
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Fig. 10. Pale Chub GSI and Korean Chub GSI during the survey period. 


ticular period. The importance of knowing about the spawn- 
ing period was indicated by Schlosser (1990), which is an 
important component in determining basic life-history infor- 
mation, for assessing the impacts of environmental variability 
on the dynamics of fish populations. 
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